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Introduction: We model the interior of Rhea 

based on observational constraints and the results from 
geodynamical models available in the literature. Ten 
main types of models are defined, depending on the 
presence or absence of a high-pressure ice layer (ice 
II), and the extent of separation of the rock component 
from the volatiles. We present degree-two gravity 
components computed for each of these models in or- 
der to assess which properties of the interior are likely 
to be inferred from Cassini radio science measure- 
ments scheduled on November 26, 2005 [1]. 

State of Knowledge: Rhea is a medium-sized icy 
satellite with a radius (R) of 7 1 8 km and a mean den- 
sity of 1240±36 kg/m 3 [2]. Photogeology observations 
by Voyager show no recent geological activity ex- 
pressed at its surface [3] and high crater density versus 
little resurfacing (especially for craters smaller than 30 
km in diameter). The presence of smooth plains mate- 
rial close to the equator indicates that the satellites has 
undergone partial resurfacing in a later stage. Com- 
pressive features are also identified by [4] as 
megaridges and megascarps developed posterior to the 
formation of the cratered terrains. 

Thermal evolution models published in the ‘80s 
conclude that the energy budget of Saturn’s medium- 
sized satellite does not allow for important melting of 
the satellite ( e.g ., [5, 6]), apart from the outermost 10 
to 100 km [5]. Models including ammonia hydrates 
and salts of chondritic origin mixed with water (e.g., 

[7, 8]) conclude to more melting and advanced partial 
differentiation. 

We consider ten different models that differ by (a) 
the extent of differentiation, (b) the rock phase density, 
ranging from 2700 kg/m 3 to 3600 kg/m 3 , depending on 
the degree of hydration, (c) the presence of absence of 
ice II. The pressure and temperature [5] inside Rhea 
are suitable for ice II to crystallize and this might be 
in part responsible for compression features identified 
at the surface of the satellite [4]. However, the rela- 
tionship between these features and ice II crystalliza- 
tion is not clearly understood, and the presence of this 
high-pressure ice component has been questioned by 
[5]. Models are presented in Figure 1. These models 
do not include a metallic core as thermal evolution 
studies concur that not enough energy was available 
for separation of metals from the rock component. 
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Fig 1. Models of Rhea’s interior. 

Interior Modeling: We compute ice I density pro- 
file after the equation of state presented in [9] using 
temperature profile calculated by [5]. As ice II equa- 
tion of state if poorly known, we set its density to 1200 
kg/m 3 , based on the different sources of information 
available in the literature. Self-compression in the rock 
phase component is negligible. Further details about 
Rhea’s interior modeling is presented in [10]. 

Results: Degree-two gravity coefficient C 22 are 
computed assuming hydrostatic equilibrium, as sug- 
gested by [1]. We plot C 22 as a function of the rocky 
core radius for models containing ice II and for a rock 
density of 3200 kg/m 3 in Figure 2. Results are com- 
pared to the predicted errors on the measurement of 
C 22 simulated by [1]: an absolute accuracy of lxl 0" 6 to 
7x1 0" 6 , i.e., a relative accuracy of 0.5 to 3.5%. Figure 2 
also encloses the error bar on the calculation of C 22 , 
due to uncertainties on the density profile. The latter 
are due to the poor knowledge of the equation of state 
of ice II and of the temperature profile. 

• If 2.63X10" 4 < C 22 < 2.80x1 IT 4 : the body is undiffer- 
entiated. This implies that there has not been enough 
energy available during the history of this satellite to 
allow for melting of the ice phase and the conse- 
quent segregation of the rocky material, and to bal- 
ance heat transfer through subsolidus convection. 
This also provides some constraint on the presence 
and the role of ammonia hydrates and other con- 
taminants. 
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• If C 22 is close to 2.80xl0" 4 : the satellite is homoge- 
neous, apart from slight density increase due to 
self-compression. This implies that conditions are 
not suitable for ice II to be stable. 

• If 1.67x10" 4 < C 22 < 1.90xl0" 4 : this is the indication 
of the presence of a rocky core, whose radius can 
be determined from the satellite’s mass and ices 
densities, for a given temperature profile. The un- 
certainty in the measurement of C 2 2 of 4xl0" 6 leads 
to an uncertainty in the determination of the core 
radius of 20 to 90 km. 

• If C 2 2 < 1.67xl0" 4 : this indicates either that the 
density of the rock phase is higher than the maxi- 
mum value considered in this study, or that there is 
separation of iron from the rock phase to form an 
iron core. 

• In the interval between the extreme values 
(1.67xl0 -4 and 2.80x1 0" 4 ) , we have a large uncer- 
tainty in the characteristics of Rhea's internal struc- 
ture, unless we choose an a priori model for the 
satellite. Further assumptions on the rock phase 
density and the presence or absence of ice II might 
be necessary to narrow the field of possible mod- 
els. 

We expect other observations obtained by Cassini 
to provide constraints on the nature of the rock com- 
ponent of the Saturnian satellites ( e.g ., [11]), al- 
though it is possible that it has evolved from its 
original state, for example through hydrothermal 
metamorphism. Quantification of past extensive and 
compressive geological processes inferred from high- 
resolution imaging by Cassini will help understand 
endogenic processes and thermal evolution, such as 
the succession of partial melting and refreezing 


events, as well as the potential role of volatiles such 
as ammonia hydrates. 

Conclusion. Constraining Rhea’s interior from 
gravity measurements is crucial for understanding the 
thermal evolution of this satellite. Based on current 
information, most of the ten models cannot be distin- 
guished from each other. However, assumptions on 
the density of the rock phase, presence or absence of 
ice II, and the degree of differentiation could allow a 
unique model to be determined in many cases. 

These new information for Rhea will be of high in- 
terest for understanding the interior and thermal evolution 
of the other medium-szed satellites of Saturn. 

This work is funded by the National Research Council. The re- 
search described in this (publication or paper) was carried out at 
the Jet Propulsion Laboratory, California Institute of Technology, 
under a contract with the National Aeronautics and Space Admini- 
stration. 

References: [1] Anderson, J.D. et al. (2003) Phys. 
Earth Planet. Int ., 136, 201-213. [2] Jacobson, R.A. (2004) 
Astron. J., 128, 492-501. [3] Morrison, D. et al. (1986) in: 
Satellites. Tucson, AZ, University of Arizona Press, 764- 
801. [4] Thomas, P.G. (1988) Icarus , 74, 554-567. [5] 

Ellsworth, K. and G. Schubert (1983) Icarus , 54, 490-510. 
[6] McKinnon, W. (1998) in: Solar System Ices , Publisher: 
Dordrecht Kluwer Academic Publishers, Astrophysics and 
space science library (ASSL) Series, 227, 525-533. [7] 
Consolmagno, G.J. (1985) Icarus , 64, 401-413. [8] Ste- 
venson, D. (1982) Nature , 298, 142-144. [9] Lupo, M.J. 
and J.S. Lewis (1979) Icarus , 40, 157-170. [10] Castillo, J., 
submitted to J. Geophys. Res. [12] Lunine, J. and L. Soder- 
blom (2002) Space Sci. Rev., 104, 191-208. 



Error bar on the 
calculation of 
due to the uncertainly 
on the density profile 


Same as above. 


The following features are plotted for models III and IV in which 
the rock phase density equals to 3200 kg/m 3 

— — — , Density of the ice-rock mixed layer 

Bottom radius of the outer icy shell for models with ice II 


Error bar on the measurement of equal to 0.5% 


Error bar on the measurement of C 22 equal to 3.5% 


Same as above. 


Fig. 2. Variations of C 2 2 as a function of the core radius and the bottom radius of the outer icy shell. 



